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A LIMITING ANALYSIS FOR EDGE EFFECTS IN ANGLE-PLY LAMINATES 1 


Peter W. Hsu A Carl T. Herakovlch 
Department of Engineering Science A Mechanics 
Virginia Polytechnic Institute A State University 
Blacksburg, Virginia 24061 


ABSTRACT 

This paper develops a zeroth-order solution for edge effects In angle-ply 
composite laminates using perturbation techniques and a limiting free body 
approach. The general method of solution for [te] laminates Is developed and 
then applied to the special case of a [±45] s graphite/epoxy laminate. Inter- 
laminar stress distributions are obtained as a function of the laminate thick- 
ness-to-wldth ratio h/b and compared to existing numerical results. 

The solution predicts stable, continuous stress distributions, determines 
finite maximum tensile Interlaminar normal stress o z for both [te] $ and [7e] $ 
laminates, and provides mathematical evidence for singular Interlaminar shear 
stresses t xz and i yz . 


Introduction: 

Recent numerical [1-6] and experimental [7-10] Investigations have 
demonstrated the free edge effect In composite laminates subjected to remote 
tension. Such effect has been suggested to play the dominant role In the de- 
lamlnatlon failure Initiation of some laminates. In an attempt to obtain more 
accurate free edge stress Intensltltes the problem of uniaxial extension of thin* 
elastic, balanced, symmetric, bidirectional laminates was Investigated In an 
earlier paper [11] based upon a perturbation analysis [12.13]. A key feature 
of the analysis was the force and moment equilibrium of a limiting free body 
containing the Interfaclal plane between two layers. The Interlaminar stresses 
thus obtained were compared with the finite difference solution of Pipes [6]. 

It was shown that the perturbation solution provides better results for the 
stress behavior near the free edge of the laminate. 

The present paper presents a similar analysis for angle-ply laminates by 
perturbing the three coupled dimensionless partial differential equations 
resulting from a displacement formulation. 

GowuUng Equatiom 

For the balanced, symmetric 2m layer laminate of Fig. 1, the displacement 
functions take the following forms [2]: 

u ■ c x x ♦ U(y,z) (a) 

v ■ V(y,z) (b) (1) 

w ■ W(y,z) (c) 

where c x Is the applied axial strain and U(y.z), V(y,z), and W(y,z) are three 
unknown functions. 

The dimensionless displacement equilibrium equations (with zero body forces) 
[14] are 


1 



{Q M (h/b) 2 U,YY ♦ 055 u »zz * ^26^ b ^*yy * ^ 4sV *ZZ 

♦ (Q36 ♦ Q45H h / b ) w »YZ )<k) * 0 

(02 6 (h/b)2u^ Y ♦ Q45 u *zz 4 * Q 44 V *2 Z (2) 

♦ (Q 23 ♦ Q44H h / b ) W »YZ ))(,t) * 0 

{ (Q45 ♦ Q 36 )(h/b)U iYZ ♦ (Q44 ♦ Q 23 >< h/b > v »Yz 

♦ Q44(h/b) 2 W,YY ♦ 033^*22 " 0 

where with c! 1 }* being the transformed stiffness coefficients 

of the meterlel properties from the neturel coordlnetes to the xy coordinates. 


and the largest stiffness coefficient of the k tb layer, r 

•4WC 

are the dimensionless coordinates, and U ■ U/h, V ■ V/h, and w * 

■ y/b, L ■ z/h 
> w/h are the 

dimensionless unknown displacement functions. 

Symmetry conditions lead to [14] 

U(Y,Z) - U(Y,-Z) 

(a) 


V(Y,Z) • V(Y,-Z) 

(b) 


W(Y,Z> • -W(Y.-Z) 

(0 

(3) 

U(Y.Z) ■ -U(-Y.Z) 

(d) 

V(Y.Z) - -V(-Y.Z) 

(•) 


W(Y,Z) - W(-Y,Z) 

(f) 



which yield the following symmetry constraints on the displacement functions: 


{U, 2 (y,0)> (m) 

• 0 

{V, 2 (y,0)) (m) 

- 0 

{W(y,0))^ m) • 

0 

{U(0,z)} (k) ■ 

0 

{V(0,z)) (k) ■ 

0 

{W, y (0,z)) (k) 

- 0 


(4) 


( 5 ) 


4 


4 


if 
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where m denotes the layer adjacent to the midplane Z • 0, and arbitrary layers 
are denoted by k. 

The appropriate stress free boundary conditions can be expressed as: 

o^ttt.Z) • ^ V. y (i1.Z) ♦ W.j(il.Z) ♦ ^ U, y (tl.Z)) <k) -0 (a) 

C C c c c 

t*y*< 4l * Z > * -r v -yC l1 ■" ♦ -jr W.jtil.Z) ♦ -§* U. y (tl.Z)) (k) .0 (b)(6) 

t* k) (*i,z) • v. z (*i.z) ♦ ^ w. y (ii.z) ♦ ^ u, 2 (ii,z)) <k) -o (c) 

along the free edges, and 

o ( J ) (y.±l) - ^ V,y(Y,±l ) ♦ ^ W, z (Y.t1) ♦ ^ u. y (y.±1)) (1) - 0 (a) 

c c c 

.<’>(», ±1) ■ <-jp *. Z (».»U ♦ -£■ W, y (V,t1) ♦ -j!£ U, z (y.il)) (,) - 0 (b)(7) 

■ < TT V, Z (y,,,) + T S (Y,l1) * ^ U. z (Y.±m (1) - 0 (c) 

on the top and the bottom surfaces. 

To solve the boundary value problem defined by Equations (2) and (4)-(7), 
only the first quadrant of the YZ-plane needs to be considered due to the favor- 
able symmetry of the laminate. Recognizing that Lhe boundary layer effect exists 
near the free edge Y ■ 1, the perturbation solution Is sought by considering two 
regions of the laminate: the Interior region (away from the free edge) and the 

boundary layer region (near the free edge). 

PeAluAbation Solution 

(1) The Interior region (0*Y<1) 

In this region the free edge stress boundary conditions (6) are dropped and 
attention Is focused on the solution to Equations (2) subject to Equations (4). 

(5) and (7). To seek a straight forward asymptotic expansion, let 


■■ t ■ ■ —<— mm .. . ■ 4 — ■ i i|i t — — I 

4 

U< k > • E c"U n (k) ( V * Z > 
n»0 

V (k) • E c n V n (k) (Y.Z) 
n«0 

W< k > - E c n W n (k) (Y.Z) 
n«0 

where tho subscript 1 denotes the Interior region and the small parameter c(«1) 
represents the th1ckness*to-w1dth ratio h/b. Substituting these expansions Into 
Equations (2) and equating coefficients of equal powers of c to zero result In 
Infinite sets of equations. The zeroth-order equations take the form 

c ° : {^ 55 u o,ZZ * ^45^0, ZZ 
|Q45 u o.ZZ + ^44 v o,ZZ 
|Q33 w o,ZZ 

As a result of the symmetry conditions (3), Equations (9) have solution 
In the form 

U 0 (k) ■ B 0 (k> (Y) 

V 0 (k) ■ D 0 <k) (Y) (10) 

W < k > ■ E (k) r 
0 0 

where B^(Y), D^(Y), E^ must satisfy the vanishing stress conditions to 
recover the lamination theory In this Interior region. That Is, 

c c c 

(C, 3 { x ( l-e)Y/h ♦ -jp D 0 (Y) ♦ -jS E fl Y ♦ -j£ B 0 (Y)) (k) • 0 


< k > .0 
< k > . 0 
(k) ■ o 


(9) 


(e) 


k • 1,2,3,... ,m 




( 11 ) 


t 


s 

which irt In effect generalizations of Equations (7a) which satisfy the symmetry 
constraints (4) and (5). The nonzero centrjl plane stresses are now required to 
satisfy the equlllbrlua conditions (Flq. 2): 

k ",(°i k)(0 - 2,h k)‘ " * ¥ E o )(k)h k { « 

♦ " C ? L ^ h k O 0>>(») ♦ * — B< o )<¥ )| “ 0 (,2) 

k-1 6 0 kO b 0 ' 

k * ) (,<5>(o.i)hk)- h|j j (e )t o.«)/h * t 0 )« k) v , 

♦ * ^i|i!k B (k>( V ) * ■ 0 (,3) 

Enforcing exact dlsplaceaent continuities In U and V across each Interface, 
yields 

B 0 (1) (Y) • B 0 (2) (Y) ■ .... ■ B 0 (k) (Y) • ff 0 (Y) (14) 

D 0 ^)(Y) • D 0 (2 '(Y) - .... ■ 0 o (k) (Y) - b' 0 (Y) 05) 

These equations reduce Equations (11)-(13) to m+2 equations for the m+2 unknowns 
B 0 (Y). 0 o (Y), E^. The solution to these reduced equations uniquely determines 
the zeroth-order Interior region solution (10). This Is the solution from 
lamination theory as will be demonstrated later for a four layer angle-ply 
laminate. 

It Is Important to note that although exact displacement continuity In W was 
not Imposed In the Interior region. It will be shown to be satisfied automatically 
for angle-ply laminates. 
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(2) The boundary layer rag Ion (Y«1) 

Introducing tha itratchlng transformation 

n ■ (l-Vl/g) 

naar tha fraa adga Y ■ 1 to tha govamlng aquations (2) rasults In tha following 
equations: 

<Q 66 U, nn + Q 55 U ’7.Z ♦ ^'nn * °45 V »ZZ " (Q 36 * 045> w 'nZ> (k) " 0 

<Q 26 U. nn + W’ZZ ♦ <>22 v »nn * W'll ” < Q 23 * ^44> w *nZ> ( k 5 ' 0 (b) ° 7) 

(k) 

{-(O 45 ♦ Q36) u » n Z * W 44 * ^23^’nZ * ^44 w »nn * ^33 W *Z2 } “ 0 

To soak a solution which satisfies tha symmetry conditions (3), tha constraint 
equations (4), and the asymptotic recova.y of tha lamination theory for large n, 
the following expansions are assumed: 

■‘S’ ’ a nm ♦ v v ~ *' 

v< b’ ' „. E 0 [0 B (t) ♦ R„« X| ' n C0$ . n z] !k) t n (18) 

M ( S ) • j g ty * v V i,n V ]<l ‘’ c " 

where D^, E^ are the Interior region solution given by Equations (11)- 
n n n 

(15), P^, and are undetermined coefficients, and </*' are undetermined 
positive quantities (In radians). The subscript b denotes the boundary layer 
region. 

Substituting Equations (18) Into Equations (17) and neglecting higher-order 
terms results In the following set of three simultaneous algebraic equations 
corresponding to the order c°: 
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^66 x o* " ^55 ®o 2,P o 4 ^26 x o ? " ^45°o ) R o " ^36 4 ^45^ x o a o s o^ ^ " 0 

C(026 x o 2 * Q 45 °o 2 ^ P o 4 ^2 o Z * ®44^# 2 ^ “ ^23 4 ‘WV’oV^ " 0 vl9) 

((Q45 ♦ 036) x o a o p o 4 ^44 4 ^23^ x o°o r o 4 ^44 x o 2 * ^33“o 2 ^o^ ^ ' 0 

k ■ 1,2 m 

For each nontrivial tarn of Equations (18) to axlst tha determinants of the 
algebraic aquations (19) must vanish Individually. This leads to a sixth-order 
algebraic equation for each layer which can be regarded as a third-order equation 
by classical treatment [15]. The six roots may be expressed as 

U 0 (1 ,2) ■ i a uo) (k) 

(X 0 (3,4) ■ i F a 0 ) (k) 

fx 0 ( 5 . 6 ) ■ t c o 0 )( k) 

where 6^ k \ c^ are three constants In terms of the material 
properties of the k th layer. For matching considerations, however, the 
positive roots must be dropped since they lead to exponential growth of the 
displacement, strain and stress fields for large n (or small Y). Hence, the 
zeroth-order expansions of Equations (18) take the following general form: 

_ _ - (It) 

U<fc> • (B 0 (Y) ♦ <P,«“ a °" ♦ ♦ P 3 '‘ Ca ° n > C0 ‘ »o I> 

y(];) . t0 o (Y) ♦ (R,e au ° n ♦ R2« f “° n ♦ R 3*''’° n ) e0 * “o z) * * 
b 0 1 

„(k) . , EoZ ♦ ( S ,i*“on * ♦ s/"®") «tn « 0 z »' k ’ 


( 21 ) 
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where P Q has been replaced by Pj, P ? , Pj, ttc. 

It nty be shown that duo to tho sops rate variable nature of Equations (21), 
no exact satisfaction of the free edge stress boundary conditions (6a) and (6b) 
(for all Z) and the stress boundary conditions (7) on the top and bottom sur- 
faces (for all Y) can be achieved. By arguing that higher-order term:, serve 
as correction terms, attention can now be focused on points (n"0,Z k +c) and 
(i^OtZ^-c) on the free edge (Fig. 3). That Is, requiring exact satisfaction 
of the boundary conditions (6a,6b,6c) at these points only [11,14] with 
Equations (21) and considering the resulting force and moment equilibrium 
of the limiting free body of thickness 2; (0<c«<1) result In the algebraic 
equations: 

{[C 2 6^ p i^ p 2^ p 3) * C 22 (aRj*6« 2 *cR 3 ) 

♦ C 23 (S 1 +S 2 ^S 3 )]o 0 cos(a 0 (Z k tc)) («) 

* < c \2^ x ♦ ir E o * t 2 - °i <‘» ♦ T* "o <‘»]*> tl0 

{£C$ 6 (* p 1 + 6 p 2 +c p 3 ) ♦ C^SRj+BR^cRj) 

* C 36 ( S 1 * V S 3 ) \ cot (« 0 < V c * ) (I b , ) ' ( : 22 ) 

• • t 'leV 1 * * T E o ♦ ^ “o' <*’> * ¥ B 0 

IC 44 [(R 1 « 2 « 3 ) - (S,J*S 2 B.SjC)] ♦ C 45 (P,+P 2 *P 3 )) (k) ■ 0 (c) 

k ■ 1 ,2 

Note that the right hand sides of Equations (22a) and (22b) are all known 
quantities from the Interior region solution. Solving nine equations (three 
from (17a), three from (17b) and three boundary conditions (22a)-(22c)) leads 
to the determination of the nine unknown coefficients Pj, Rj, (1*1,3) In 


term of a^. The validity of the solution thus obtained can be readily 
checked by the self-equilibrating condition of the stress resultant 


/: 




0. 7 dY r 0 


( 23 ) 


for any level of Z. Finally, equating the force and moment resulting from the 
boundary layer displacement fields (21) with the Interior region stress 
resultants determines the values of o^((Z k tc)) to their order of accuracy. 
Thus, the "near* Interlaminar" stress distributions can be obtained based upon 
a reference layer. It should be noted that displacement continuity In this 
boundary layer region has not yet been Imposed as a physical requirement. 

It will be Imposed subsequently In a numerical example. 

Foua Laytx Angli-Ply Laminate a 

For advanced fiber-reinforced composites having three mutually perpendlcu- 
lar planes of elastic symmetry, the stiffness coefficients C^ 7 vanish. For a 
four layer angle-ply laminate with synvnetrlc [±e] orientations (Fig. 4a), the 
following relations between material constants (with respect to xyz coordinates) 
are found to exist [14] 


c u” 





1 ■ 1 ,2,3 and J ■ 1 ,2,3 
k - 4,5,6 
n ■ 1 ,2,3 


The zeroth-order Interior region solution (10) yields 
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■ «<*> ■ 0 

0 Ml 

. , ») . . . ( . c ,» c » ' Wn' 

o < c i2 c 33 - c 23 c 23) 1 *** 



(2) a _ < c 13 c 22 


( C 22 C 33 


C 12 C 23^ ] 

TTT 


Id-c) 


C 23 C 23 




(24) 


It can be teen from (24) that the exact continuity In W results automatically. 

On the central plane (Y«0), the stresses are obtained by combining 
Equations (24), the constitutive equations, and the strain-displacement re- 
lations. The results are 

Oy (,, (0.Z) * Oy (2) <0,Z) • 0 (,) 


T X y (1 >(0,Z) ■ - T xy < 2 >(0,Z) 


(25) 


[ c . CZ6< C 1Z C 33 - C 13 C 23> * C 36 (C 13 C 22 ' (1.,)(h) 

L 6 ( c 22 c 33 * c 23 c 23^ J 


which Indicates that the lamination theory (or the zeroth-order interior 
region solution) contributes no normal stress along the central plane Y-0. 
For equilibrium considerations the Interlaminar shear stress resultant and 
tV couple moment due to the Interlaminar normal stress a z should both be 
expected to vanish (Fig. 5). Thus, two more self-equilibrating conditions 
are established. In addition to Equation (23), as 


Jo v * ■ 0 
f a z ydy • 0 


(26) 


( 27 ) 
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Mhere t y2 and a z are both determined by solving tht boundary layer equations 
(16-22) with k-1,2. 

Numerical exsmplo 

(1) [45/-45] # graphite-epoxy laminate 

Consider the [45/-45] # graphite-epoxy laminate of constant layer thickness 
h/2 (Fig. 4a). The material properties are 


E n - 20 x 10° (psl) 

E 22 " E 33 * 21 x 1 ° 6 <P* 1 ) 
G 12 * G 23 " 6 13 “ 0,85 * 

v 12 " v 23 “ v 13 " 0,21 

The transformed stiffness coefficients are 


4S(x 10~ 6 psl) 


c„ (,) 

■ 

6.745 

c (2 > 
c n 


m 

5.045 

c (2) 
C 12 

c (,) 

c 13 

m 

0.521 

c‘ (2) 

C 13 


0) 
72 

(D 

'23 

(D 

'33 

( 1 ) 

'16 

(1) 
'36 
(1) 

'44 

, 0) 

“56 

c 0) 

L 45 


C 2 

C 2 

C, 

% 

c , 

C, 

% 

C* 


6.745 

0.521 

2.213 

■ -4.506 


-0.04387 


C 0) 
c 55 

5.33 

0 


0 So 


( 2 ) 

22 

( 2 ) 
23 
( 2 ) 
33 

( 2 ) 

16 

. ( 2 ) 

'36 

. ( 2 ) 

'44 


C, 

C 4b 


( 2 ) 


“66 

r (2) 


•45 (x 10“ 6 p>1) 
6.745 

5.045 

0.521 

6.745 

0.521 

2.213 

C 2 ^ 2) - 4.506 
0.04387 

C 5 r ■ °-85 

5.33 
0 


(28) 


(29) 


12 


The Interior solution (24) yields 


u o (1) B U 0 (2) " 0 

V o O) . v 0 C« ■ -0.7433 c x Jf V(l-e) 
Mo (1) ■ Ho ( 2) • -0.0604 c x Z (1-c) 

which lead to the central plane stresses 

T il > (0»Z) - - t< 2) (O.Z) - 1.154 cjl-c)(10 6 ps1) 


( 30 ) 


xy 


xy 


o^’^O.Z) • - o< J) (0,Z) • 0 


(31) 


The boundary layer equation (17-19) yield two Identical sixth-order algebraic 
equations [14] which give three pairs of real roots. For matching considerations 
only the three negative roots are taken. Finally the composite solution (In 
the perturbation sense) Is formed as 

U c (lc) ■ ((P^* 1 0 ♦ Pf 2 ° n ♦ P 3 e* 3 ° n )cos a Q Z}^ 

V c (k) - - 0.7433 C x (l-c)£ T ♦ * R 2 e 02a ° 

♦ Rj^^Jcos oi 0 Z) (k) 

« c < k >-. 0.0604 C x (1.«)2 M(V V 2V 

♦ S 3 e 63< *° n )s1n a 0 Z) (k) 


(32) 


where 


*1 


00 

00 


1.2364 

0.2903 


6 3 


(k) . 


0.9659 
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Exact satisfaction of the governing aquations and tha boundary conditions at 
points (n*0»1/2+c) and (n-0,1/2-c) loads to tha following aquations: 


p/ 1 ) - - 0.5871 
P 2 (1 > - 0.1707 a 1 

PjO) ■ 1.2021 * } 

r/ 1 ) - - 0.6309 *1 
• - 0.1813 ^ 
R 3 ^) ■ 1.1897 

S,0) - 1.1358 ♦ 1 

S-^ 1 * ■ 0.0347 S] 

S 3 O) . - 1.0736 ^ 


, p/ 2 > ■ 0.5871 + 2 

P 2 (2) - - 0.1707 s 2 
P 3 (2) - - 1.2021 4 2 
R 1 (2) - - 0.6309 « 2 
R 2 (2) - - 0.1813 * 2 
R 3 ( 2 ) - 1.1897 + 2 

S / 2) ■ 1.1358 * 2 

S 2 < 2 > ■ 0.0347 + 2 

, S 3 (2) * - 1.0736 t 2 


(33) 


where 

* x (l-e) 

a 0 tl, COS(a 0 (1, (l + ;) ) 

* 2 a 0 (Z, COS(a 0 lZ, (J. - 0) 

0 < c «< 1 


(34) 


It can be shown that these coefficients lead to Identical satisfaction of the 
self-equilibrating conditions (Equations (23) and (26)) when the lower limit 
Is replaced by Infinity - the corresponding zeroth-order domain of the 
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Inttrlor region. Furthermore, requiring the force equilibrium condition 
f t<‘> (y. 1/2 t 0 <ty ♦ t<J> (0,i) $ • o 

ad 'i l i u M ) . 


leads to 


(k) 


0.5 (0 < c «< D 


( 35 ) 


Now consider Equations (32) and (33). It Is clear that layer 1 (>45° ) 
and Layer 2 (-45°) are antisymmetric In U and symmetric In both V and W with 
respect to the Infinitesimal thin slice (Fig. 4d). Upon enforcing exact 
continuity In displacement U, V, W at l ■ 1/2, the following equation Is 
obtained 


11m cos 
c-*o 


(D 


♦ a 0 0) C 


11m cos 

C-K) 


(2) 


4 - - 


which yields 


f 0) 


cos 


♦ c . <’> 


( ( 2 ) 


■ cos 


2 c “o 

for 0 < c «< 1 


( 2 ) 


(36) 


(37) 


and 


11m 

C-K> 




( 38 ) 
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The Inter laminar s trasses based upon the lower layer (-45°) may now be ex- 
pressed from the stress-displacement relations as 

t x2 ■ (0.85 x 10 6 )(l-c)t x [-0.5871i 1,2364an ♦ 0.1707e 0,2903aO 

♦ 1.2021e°- 9659an ]tan(| - ac) (a) 

t _ - (0.85 x 10 6 )(l-e)t y [2.0350a 1 ,2364an ♦ 0.1913e 0,2903an 

- 2.2263e 0,9659on ]tan(£ - ac) (b) (39) 

°z " O-Ot x (10 6 )[2. 1389a 1 * 2364an ♦ 0.0472e 0,2903an 

- 1.8281e 0,9659an ] (c) 

0<c«<l 

If the stacking sequence Is reversed to [-45/45] $ , (Fig. 4e) the Interlaminar 
stresses become 

t xz • (0.85 x 10 6 )(1-c)c x [+0.5871e 1,2364an - O.^e 0,2903 ® 0 

- 1.2021e 0,9659an ]tan(§- - ac) (a) 

r - (0.85 x 10 6 )(1-c) t [2.0350e 1,2364an ♦ 0.1913e 0,2903an 
yz a 

- 2.2263e 0,9659an ]tan(| - ac) (b) (40) 

°z " 0-Oc x (10 6 )[2. 1389c 1 ,2364an + 0.0472e°* 2903ari 

- 1.8281e 0,9659on ] (c) 

0<c <<< 1 

Results and Discussion 

From Equations (39) and (4G) It Is clear that the Interlaminar shear 
stresses t x2 and t yz are both proportional to the near singular value of 
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tan(J - ac) which results frost Equation (35 6 38). Hone* a serves as a problem 
parameter which may be more realistically determined experimentally. Figure 
6 shows the Influence of a on the Interlaminar shear stress t x2> Obviously. 
t X 2 becomes more singular as a Is Increased and may attain a much higher 
finite maximum value at the free edge than the calculated finite difference 
result of [3]. This Is In agreement with the work of Pipes and Pagano [7] 

In which they found that x %z tends to grow without bound. Figure 7 shows 
the Interlaminar shear stress T yz as a function of the problem parameter 
a. Although Ty Z Is proportional to the near singular value tan(| -a;) 

(0<c«<l) It Is zero at the free edge thus satisfying the stress free boundary 
condition. It should be noted that i attains larger peak values, for higher 
values of o. The negative-positive variation of T yz confirms the validity of 
the self-equilibrating condition (26). The finite difference solution, however, 
does not predict such variation. In a later paper, It will be shown that the 
negative-positive variation agrees well with the finite element result by 
Renlerl [16] which further supports the present theory. The variation of the 
Interlaminar normal stress o z In Fig. 8 (Eon. 40c) Indicates that the maximum 
finite value of o z at the free edge Is Independent of the problem parameter a. 
The only Influence of a on o z lies In the boundary layer width. The positive- 
negative (tensile-compressive) variation of o z confirms the solution validity 
by satisfying the self-equilibrating condition (23). The finite difference 
results, on the other hand. Indicate Instability near the free edge [14]. 

The present theory (Eqns. (8) - (22)) Is based upon the zeroth-order 
analysis of the geometric ratio h/b. Hence the smaller h/b, the better the 

* All finite difference reuslts presented In this paper were obtained by 
the authors using the program supplied by Professor Pipes. 


■ ■ ■ 1 1 1 1 1 
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solution accuracy [14]. Tha affacts of this ratio on tha Interlaminar strass 
components can ba obsarvad In Figs. 9-11. In Fig. 9. for a smaller h/b, ? 
has a smaller boundary layer width while attaining a higher maximum value at 
tha free edge (for a fixed a). Similar behavior Is found for r yz (Fig. 10) 
where tha strass attains a higher peak value and a smaller boundary layer 
width for a smaller h/b. In Fig. 11, a higher a _ and a smaller boundary 
layer width are obtained for a smaller h/b. The Interlaminar stress distri- 
butions for the reversed stacking sequence [-45/45] s Gr/E are not plotted. 
However It Is Important to point out [14] that only t X 2 experiences a sign 
reversal when the stacking sequence Is reversed. The sign of the remaining 
two components of Interlaminar stress Is not a function of the stacking 
sequence. Thus for reliable design fo angle-ply laminates, the delamlnatlon 
failure mode due to the tensile o z at the free edge should always be taken 
Into consideration. 

Conclusions 

A method of solution for the problem of elastic, balanced, symmetric 
laminates subject to uniaxial extension has been developed based upon the 
perturbation theory. Attention has been focused on the force and moment 
equilibrium for an Infinitesimally thin slice containing the Interfaclal 
plane. The solution provides better Insight Into the free edge Interlaminar 
stress behavior for thin angle-ply laminates (h/b«1) than existing numerical 
solutions. 
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FIGURE 2 FREE BODY DIAGRAM OF QUARTER YZ- PLANE 
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FIGURE 4. FOUR LAYER ANGLE - PLY LAMINATE 
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FIGURE 6. PRESENT THEORY AND FINITE 
DIFFERENCE RESULTS FOR r XI IN [45 /-45] 
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FIGURE 7. PRESENT THEORY AND FINITE DIFFERENCE 
RESULTS FOR T in [45/- 45], 
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JURE 9. INTERLAMINAR SHEAR STRESS 
FUNCTION OF h/b IN [45Z-45], 




FIGURE 10. INTERLAMINAR SHEAR STRESS 




FIGURE II. INTERLAMINAR NORMAL STRESS <r z 
FUNCTION OF h/b IN [45/-451 





